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ABSTRACT 

We study the waiting time distributions of solar flares observed in hard X-rays with ISEE- 
3/ICE, HXRBS/SMM, WATCH/GRANAT, BATSE/CGRO, and RHESSI. Although discordant 
results and interpretations have been published earlier, based on relatively small ranges (< 2 
decades) of waiting times, we find that all observed distributions, spanning over 6 decades of 
waiting times (At rj 1CP 3 — 10 3 hrs), can be reconciled with a single distribution function, 
N(At) oc Ao(l + AoAi) -2 , which has a powerlaw slope of p « 2.0 at large waiting times 
(At « 1 — 1000 hrs) and flattens out at short waiting times At < Aio = 1/Ao- Wc find a 
consistent breakpoint at At = 1/A = 0.80 ± 0.14 hours from the WATCH, HXRBS, BATSE, 
and RHESSI data. The distribution of waiting times is invariant for sampling with different 
flux thresholds, while the mean waiting time scales reciprocically with the number of detected 
events, Ato oc 1/ndet- This waiting time distribution can be modeled with a nonstationary Pois- 
son process with a flare rate A = 1/At that varies as /(A) oc A -1 exp — (A/Ao). This flare rate 
distribution represents a highly intermittent flaring productivity in short clusters with high flare 
rates, separated by quiescent intervals with very low flare rates. 

Subject headings: methods: statistical — instabilities 

1. INTRODUCTION 

You might drive a car in a foreign city and have to stop at many red traffic lights. From the statistics 
of waiting times you probably can quickly figure out which signals operate independently and which ones 
operate the smart way with inductive-loop traffic detectors in a so-called demand- actuated mode. Thus, 
statistics of waiting times bears crucial information how a system works, either having independent elements 
that act randomly, or consisting of elements with long-range connections that enable coupling and synchro- 
nization. In geophysics, aftershocks have been found to exhibit different waiting time statistics (Omori's 
law) than the main shocks of earthquakes (e.g., Omori 1895; Bak et al. 2002; Saichev and Sornette 2006). In 
magnetospheric physics, waiting time statistics is used to identify Poisson random processes, self-organized 
criticality, intermittent turbulence, finite system size effects, or clusterization, such as in auroral emission 
(Chapman et al. 1998, 2001), the auroral electron jet (AE) index (Lepreti et al. 2004; Boffetta et al. 1999), 
or in substorms at the Earth's magnetotail (Borovsky et al. 1993; Freeman and Morley 2004). Waiting time 
statistics is studied intensely in solar physics, where most flares are found to be produced by a Poissonian 
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random process, but there are also so-called sympathetic Bares that have a causal connection or trigger each 
other (e.g., Simnett 1974; Gergely and Erickson 1975; Fritzova-Svestkova et al. 1976; Pearce and Harrison 
1990; Bumba and Klvana 1993; Biesecker and Thompson 2000; Moon et al. 2002). Waiting time statistics 
of solar flares was studied in hard X-rays (Pearce et al. 1993; Biesecker 1994; Crosby 1996; Wheatland et 
al. 1998; Wheatland and Eddy 1998), in soft X-rays (Wheatland 2000a; Boffetta et al. 1999; Lepreti et 
al. 2001; Wheatland 2001; Wheatland and Litvinenko 2001; Wheatland 2006; Wheatland and Craig 2006; 
Moon et al. 2001; Grigolini et al. 2002), for coronal mass ejections (CMEs) (Wheatland 2003; Yeh et al. 2005; 
Moon et al. 2003), for solar radio bursts (Eastwood et al. 2010), and for the solar wind (Veltri 1999; Podesta 
et al. 2006a,b, 2007, Hnat et al. 2007; Freeman et al. 2000; Chou 2001; Watkins et al. 2001a,b, 2002; Gabriel 
and Patrick 2003; Bristow 2008; Greco et al. 2009a, b). In astrophysics, waiting time distributions have been 
studied for flare stars (Arzner and Guedel 2004) as well as for black-hole candidates, such as Cygnus X-l 
(Negoro et al. 1995). An extensive review on waiting time statistics can be found in chapter 5 of Aschwandcn 
(2010). 

In this study we focus on waiting time distributions of solar flares detected in hard X-rays. The most 
comprehensive sampling of solar flare waiting times was gathered in soft X-rays so far, using a 25-year 
catalog of GOES flares (Wheatland 2000a; Boffetta et al. 1999; Lepreti et al. 2001), but three different 
interpretations were proposed, using the very same data: (i) a not-stationary (time-dependent) Poisson 
process (Wheatland 2000a), (ii) a shell- model of MHD turbulence (Boffetta et al. 1999), or (hi) a Levy 
flight model of self-similar processes with some memory (Lepreti et al. 2001). All three interpretations can 
produce a powerlaw-like distribution of waiting times. On the other side, self-organized criticality models 
(Bak et al. 1987, 1988; Lu and Hamilton 1991; Charbonneau et al. 2000) predict a Poissonian random 
process, which has an exponential distribution of waiting times for a stationary (constant) flare rate, but 
can produce powerlaw-like waiting time distributions with a slope of p < 3 for nonstationary variations of 
the flare rate (Wheatland and Litvinenko 2002). Therefore, the finding of a powerlaw-like distribution of 
waiting times of solar flares has ambiguous interpretations. The situation in solar flare hard X-rays is very 
discordant: Pearce et al. (1993) and Crosby (1994) report powerlaw distributions of waiting times with a 
very flat slope of p « 0.75, Biesecker (1994) reports a near-exponential distribution after correcting for orbit 
effects, while Wheatland et al. (1998) finds a double-hump distribution with an overabundance of short 
waiting times (At w 10 s - 10 min) compared with longer waiting times (At s; 10- 1000 min), but is not 
able to reproduce the observations with a nonstationary Poisson process. In this study we analyze flare 
catalogs from HXRBS/SMM, BATSE/CGRO and RHESSI and are able to model all observed hard X-ray 
waiting time distributions with a unified model in terms of a nonstationary Poisson process in the limit of 
high intermittency. We resolve also the discrepancy between exponential and powerlaw-like waiting time 
distributions, in terms of selected fitting ranges. 



2. Bayesian Waiting Time Statistics 

The waiting time distribution P(At) for a Poissonian random process can be approximated with an 
exponential distribution, 

P(At) = Ae- AAt , (1) 

where A is the mean event occurrence rate, with the probability / P(At)dAt = 1 normalized to unity. If the 
average event rate A is time-independent, we call it a stationary Poisson process. If the average rate varies 
with time, the waiting time distribution reflects a superposition of multiple exponential distributions with 
different e-folding time scales, and may even ressemble a powerlaw-like distribution. The statistics of such 
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inhomogeneous or nonstationary Poisson processes can be characterized with Bayesian statistics. 



A nonstationary Poisson process may be subdivided into time intervals where the occurrence rate is 
constant, within the fluctuations expected from Poisson statistics, so it consists of piecewise stationary 
processes, e.g., 

A ie - AlAt 
A 2 e-^ At 



P(At) 



for t 1 <t<t 2 
for t 2 <t<t 3 



(2) 



A, 



-A„ At f or t < £ < £ 



n+1 



where the occurrence rate Aj is stationary during a time interval [t i7 t i+ i], but has different values in subse- 
quent time intervals. The time intervals [U,ti + i] where the occurrence rate is stationary are called Bayesian 
blocks, and we talk about Bayesian statistics (e.g., Scargle 1998). The variation of the occurrence rates 
Ai, A2, A„ can be defined with a new time-dependent function X(t) and the probability function of waiting 
times becomes itself a function of time (e.g., Cox & Isham 1980; Wheatland et al. 1998), 



P(t, At) = X(t + At) exp 



t+At 



X(t')dt' 



(3) 



If observations of a nonstationary Poisson process are made for the time interval [0, T], then the distribution 
of waiting times for that time interval will be, weighted by the number of events X(t)dt in each time interval 
{t,t + dt), 



P(Ai) = l J X(t) P(t,At)dt, 



(4) 



where the rate is zero after the time interval t > T, i.e., X(t > T) = 0, and N = J Q X(t)dt. If the rate 
is slowly varying, so that it can be subdivided into piecewise stationary Poisson processes (into Bayesian 
blocks), then the distribution of waiting times will be 



P(At) «^(A 2 )A,e- 



A;At 



(5) 



where 



</>(Ai) 



XAj 



(6) 



is the fraction of events associated with a given rate A^ in the (piecewise) time interval U (or Bayesian block) 
over which the constant rate Ai is observed. If we make the transition from the summation over discrete time 
intervals U (Eqs. 5 and 6) to a continuous integral function over the time interval [0 < t < T], we obtain, 



/A , C X(t) 2 e~ x ^ At dt 

P(At) = ^ — y -4 

fix(t)dt 



(7) 



When the occurrence rate X(t) is not a simple function, the integral Eq. (7) becomes untractable, 
in which case it is more suitable to substitute the integration variable t with the variable A. Defining 
/(A) = (l/T)dt(X)/dX as the fraction of time that the flaring rate is in the range (A, X+dX), or f(X)dX =dt/T, 
we can express Eq. (7) as an integral of the variable A, 



P(At) 



2„-AAt 



dX 



/ °°A/(A) dX 



(8) 
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where the denominator Ao = J °° Xf(X)dX is the mean rate of flaring. 

Let us make some examples. In Fig. 1 we show five cases: (1) a stationary Poisson process with a constant 
rate Ao; (2) a two-step process with two different occurrence rates Ai and A2; (3) a nonstationary Poisson 
process with piece-wise linearly increasing occurrence rates X(t) = Xoit/T), varying like a triangular function 
for each cycle, (4) piece-wise exponential functions, and (5) piece-wise exponential function steepened by a 
reciprocal factor. For each case we show the time-dependent occurrence rate X(t) and the resulting probability 
distribution P(At) of events. We see that a stationary Poisson process produces an exponential waiting time 
distribution, while nonstationary Poisson processes with a discrete number of occurrence rates A^ produce a 
superposition of exponential distributions, and continuous occurrence rate functions X(t) generate powerlaw- 
like waiting time distributions at the upper end. 

We can calculate the analytical functions for the waiting time distributions for these five cases. The 
first case is simply an exponential function as given in Eq. (1) because of the constant rate X(t) = A , 

P(At) = A e- AoAt . (9) 

The second case follows from Eq. (5) and yields 

P(At) = ^A ie - AlA * + ^A 2 e- A2A * . (10) 

The third case can be integrated with Eq. (7). The time-dependent flare rate grows linearly with time 
to a maximum rate of X m = 2Ao over a time interval T, with a mean rate of A , 

A(i)=A ro (|)=2Ao(!) (11) 

Defining a constant a = -X m (At/T) the integral of Eq. 7 yields P(At) = (2A m /T 3 ) / Q T ' t 2 e at dt. The integral 
J x 2 e ax dx = e ax (x 2 /a — 2x/a 2 + 2/a 3 ) can be obtained from an integral table. The analytical function of 
the waiting time distribution for a linearly increasing occurrence rate is then 



P(At) = 2X r , 



e -AoAt 



1 



(A m Ai) 3 \(X m At) (X m At) 2 (A m Ai) 3 



(12) 



which is a flat distribution for small waiting times and approaches a powerlaw function with a slope of p = 3 
at large waiting times, i.e., P(At) oc Ai~ 3 (Fig. 1, third case). The distribution is the same for a single 
linear ramp or a cyclic triangular variation, because the total time spent at each rate [A, X + dX] is the same. 

The fourth case, which mimics the solar cycle, has an exponentially growing (or decaying) occurrence 
rate, i.e., 

1 \ fx 



fW = [-J^{--J- (13 » 

defined in the range of [0 < A < 00] , and has a mean of A . The waiting time distribution can therefore be 
written with Eq. (8) as 

P(At) = j™^ cxp(-A[i + Ao At])dA, (14) 

which corresponds to the integral x 2 e ax dx = —2/a 3 using a = —(1 + X At)/Xo and thus has the solution 
P(At) = -2/(o s A5), i.e., 
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For very large waiting times (At ~> 1/Ao) it approaches the powerlaw limit P(At) « (2/Ajj)Af 3 (see Fig. 1 
fourth panel). 

The fifth case has an exponentially growing occurrence rate, multiplied with a reciprocal factor, i.e., 

/(A) = A- 1 cxp (- A) , (16) 



A 



and fulfills the normalization J °° A/(A)dA = A . The waiting time distribution can then be written with 
Eq. (8) as 

p(At) = r(^) (- j- 1 1 + a ° a *]) ^ > 



o 



which, with defining a = —(1 + AoAt)/Ao, corresponds to the integral J xe ax dx — (e ax /a 2 )(ax — 1) and has 
the limit J °° xe ax dx = 1/a 2 , yielding the solution P(At) = l/(a 2 Ao), i.e., 

P(Ai) " <ITW ' <18) 

For very large waiting times (At 3> 1/Ao) it approaches the powerlaw limit P(At) w (l/A§)At -2 (see Fig. 1 
bottom panel). 

Thus we learn from the last four examples that most continuously changing occurrence rates produce 
powerlaw-like waiting time distributions with slopes of p < 2, ...,3 at large waiting times, despite of the 
intrinsic exponential distribution that is characteristic to stationary Poisson processes. If the variability of 
the flare rate is gradual (third and fourth case in Fig. 1), the powerlaw slope of the waiting time distribution 
is close to p < 3. However, if the variability of the flare rate shows spikes like (^-functions (Fig. 1, bottom), 
which is highly intermittent with short clusters of flares, the distribution of waiting times has a slope closer 
to p w 2. This phenomenon is also called clusterization and has analogs in earthquake statistics, where 
aftershocks appear in clusters after a main shock (Omori's law). Thus the powerlaw slope of waiting times 
contains essential information whether the flare rate varies gradually or in form of intermittent clusters. 



3. DATA ANALYSIS 

We present an analysis of the waiting time distribution of solar flares observed in hard X-rays, using 
flare catalogs obtained from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI), the Compton 
Gamma Ray Observatory (CGRO), and the Solar Maximum Mission (SMM), and model also previously 
published data sets from SMM (Pearce et al. 1993), GRAN AT (Crosby 1996), and ISEE-3 (Wheatland 
1998). 



3.1. RHESSI Waiting Time Analysis 

RHESSI (Lin et al. 2002) was launched on February 5, 2000 and still operates at the time of writing, 
for a continuous period of 8 years. The circular spacecraft orbit has a mean period of 96 min (1.6 hrs), a 
mean altitude of 600 km, and an inclination of 38° with respect to the equator, so the observing mode is 
interrupted by gaps corresponding to the spacecraft night (with a duration of about 35 percent of the orbit), 
as well as some other data gaps when flying through the South Atlantic Anomaly (SAA). These data gaps 
introduce some systematic clustering of waiting times around an orbital period (sa 1.6 hrs). 
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We are using the official RHESSI flare list (which can be downloaded from the RHESSI webpage 
http://hesperia.gsfc.nasa.gov/rhessidatacenter/). We are using data from the first six years of the mis- 
sion (from February 12, 2002, to December 31, 2007), when a uniform threshold for flare detection was 
applied, while a lower threshold was applied after this period. We include all flare events that have a solar 
origin, in the sense that they could be imaged in the energy range of 12 to 25 keV. The complete event 
catalog (sampled up to Feb 1, 2010) contains n = 52,014 events, of which 12,379 have been confirmed to 
be solar. The number of confirmed solar flare events we are using from the first six years includes 11,594 
events, which corresponds to a mean flare rate of < A >« 5.5 events per day during 2002-2007. 

A time series of the daily RHESSI flare rate is shown in Fig. 2 (top panel). The mean annual flare rate 
clearly drops systematically from 2002 to 2008, but the daily fluctuations are larger than the slowly- varying 
mean trend. We calculate the waiting times simply from the time difference between two subsequent flare 
events, 

Mi =U- U-i ■ (19) 

If we plot the waiting time distribution on a log-log scale (Fig. 2 bottom right), we see that the waiting time 
distribution N(At) = nP(At) can approximately be fitted with a powerlaw function, i.e., N(At) oc At~ p , 
with a powerlaw slope of p ~ 2.0. However, there is some clustering of events above one orbital period of 
At > t or i,it ~ 1.6 hrs, as well as at the second to fourth harmonics of the orbital period (Fig. 1, bottom 
left), which are caused for several reasons. A first reason is that events that start at spacecraft night cannot 
be detected until the spacecraft comes out of the night part of the orbit, causing some clustering after one 
full orbital period. A second reason is that large events that extend over more than one full orbital period 
are double-counted in each subsequent orbit. These instrumental biases have been modeled with Monte- 
Carlo simulations in previous waiting time studies (e.g., Biesecker 1994), but instead of correcting for these 
erroneous waiting times, we will just exclude the range of At ss (0.5 — 2) x t or bi t («0.4-2.4 hrs) in the fits of 
waiting time distributions. 

Thus our first result is that the waiting time distribution of RHESSI flares is approximately consistent 
with a powerlaw distribution in the time range of At wl - 1000 hrs, with a powerlaw slope of p « 2.0. For 
a nonstationary Poisson process we expect powerlaw slopes in the range of p ss 2.0, ...,3.0 (Fig. 1), so the 
measured distribution is close to what the theory predicts for very intermittently varying flare rates (Fig. 1, 
bottom). The high degree of intermittency is indeed clearly recognizable in the data (Fig. 2, top panel), 
where the short-term fluctuations are much faster than the long-term trend. One might subdivide the flare 
activity into two states, i.e., "quiescent" and "Bare-active" states, similar to the quiescent periods (soft state) 
and active periods (hard state) of pulses from accretion disk objects and black- hole candidates (e.g., Cygnus 
X-l; Mineshige et al. 1994). We indicate quiescent periods with waiting times At > 5 hrs in Fig. 2 (top 
panel), where it can be seen that they occur in every month through all 8 years, but more frequently during 
the (extended) solar minimum (2006-2008). Thus, flare-active periods are very intermittent and not lasting 
contiguously over extended periods of time. The situation is also similar to earthquakes, where aftershocks 
appear clustered in time intervals after larger main shocks (Omori's law). 

In a next step we investigate the effect of flux thresholds in the event definition on the distribution 
of waiting times, an issue that has been raised in previous studies (e.g., Buchlin et al. 2005; Hamon et 
al. 2002). Hamon et al. (2002) finds for the Olami-Feder-Christensen model (Olami et al. 1992), which is 
a cellular automaton model for systems with self-organized criticality, that event definitions without any 
threshold lead to stretched exponential waiting time distributions, while threshold-selected events produce 
an excess of longer waiting times. We investigate this problem simply by applying various flux thresholds 
to the RHESSI flare catalog, e.g., P = 10, 100,300, 1000,3000 cts s _1 , and by re-sampling the waiting times 
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for events above these flux thresholds. The corresponding six waiting time distributions are shown in Fig. 3, 
which contain n = 11,594 events detected without a threshold, and ut = 9596, 2038, 781, 271, and 108 
events for the thrcsholded subsets. The waiting time distributions clearly show an increasing excess of longer 
waiting times with progressively higher thresholds, which is expected due to the filtering out of shorter time 
intervals between flares with weaker fluxes below the threshold. Based on the reduction in the number n 
of events as a function of the flux threshold, we can make a prediction how the mean waiting time interval 
increases with increasing threshold, namely a reciprocal relationship, since the total duration T of waiting 
times is constant, 

n 

T = Atj =< n >< At >=< tit >< At >t . (20) 

i 

Thus, from the number of events nr detected above a selected threshold we can predict the mean waiting 
time, 

< At > T = <U> < At > . (21) 

< TIT > 

Using the full set of data with < n >= 11,594 events and a mean waiting time of < At >= 0.71 hrs 
(Fig. 3, top left), we can predict the distributions and average waiting times for the thresholded subsets, 
based on their detected numbers nr using Eq. (21): < At >t= 0.9,4.0,10.5,30.4,76.3 hrs. We fit our 
theoretical model of the waiting time distribution (Eq. 18) of a nonstationary Poisson process and predict 
the distributions for the threshold datasets, using the scaling of the predicted average waiting times < At >t- 
The predicted distributions (thick curves in Fig. 3) match the observed distributions of thresholded waiting 
times (histograms with error bars in Fig. 3) quite accurately, and thus demonstrates how the waiting time 
distribution changes in a predictable way when flux thresholds are used in the event selection. 

Regarding the mean flare waiting time we have to distinguish between the theoretical model value 
Ato and the mean detected interval < At > obs . The theoretical value is calculated based on a complete 
distribution in the range of At = [0, oo]. For RHESSI data we obtain a value of Ato = 0.71 hrs. The 
observational value, which is about the total observing time span (« 5.8 yrs) multiplied with the spacecraft 
duty cycle, which is about q — 0.65 for RHESSI (based on spacecraft nights with a length of 32-38 minutes), 
and divided by the number of observed events n. So, we obtain a mean time interval of < At > o6s = Tq/n = 
5.8 x 365 x 24 x 0.65/11, 594 rs 3.0 hrs. This value is about a factor of 4 longer than the theoretical value. 
This discrepancy results from either missing short waiting times predicted by the model (since the observed 
distribution falls off at waiting times of At < 0.02 hrs, i.e., w 1 minute), or because the model overpredicts 
the maximum flare rate, and thus needs to be limited at the maximum flare rate or lower cutoff of waiting 
times. Neverthelss, whatever lower cutoff of observed waiting times or what flux threshold is used, the 
theoretical value Ato always indicates where the breaking point is in the waiting time distribution, between 
the powerlaw part and the rollover at short waiting times. 



3.2. BATSE/CGRO Waiting Time Analysis 

The Compton Gamma Ray Observatory (CGRO) was launched on April 5, 1991, and de-orbited on 
June 4, 2000. The CGRO spacecraft had an orbit with a period of t or ut ~ 1.5 hrs also, and thus is subject to 
similar data gaps as we discussed for RHESSI, which causes a peak in the waiting time distribution around 
At « (0.5 - l.l)t orUt . 

We use the solar flare catalog from the Burst and Transient Source Experiment (BATSE) onboard 
CGRO, which is accessible at the NASA/GSFC homepage |Mtp:/ / umbra.nascom.nasa.gov /batse/\ Here we 
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use a subset of the flare catalog obtained during the maximum of the solar cycle between April 19, 1991 
and November 12, 1993, containing 4113 events during a span of 1.75 years, yielding a mean event rate of 
< A >s=y 6.4 events per day. BATSE has 8 detector modules, each one consisting of an uncollimated, shielded 
Nal scintillation crystal with an area of 2025 cm 2 , sensitive in the energy range of 25 keV-1.9 MeV (Fishman 
et al. 1989). 

The waiting time distribution obtained from BATSE/CGRO is shown in Fig. 4 (middle right), which 
ranges from At sa 0.01 hrs up to ~ 200 hrs. We fit the same theoretical model of the waiting time distribution 
(Eq. 18) as for RHESSI, and the powerlaw slope of p — 2 fits equally well. Since the obtained mean flare 
rate (Ai = 0.92 hrs) is similar to RHESSI (At = 0.71 hrs), the thresholds for selected flare events seems 
to be compatible. 

3.3. HXRBS/SMM Waiting Time Analysis 

While RHESSI observed during the solar cycles #23+24, CGRO observed the previous cycles #22+23, 
and the Solar Maximum Mission (SMM) the previous ones #21+22, none of them overlapping with each 
other. SMM was lauched on February 14, 1980 and lasted until orbit decay on December 2, 1989. The 
Hard X-Ray Burst Spectrometer (HXRBS) (Orwig et al. 1980) is sensitive in the range of 20-260 kcV and 
has a detector area of 71 cm -2 . The orbit of SMM had initially a height of 524 km and an inclination of 
18.6°, which causes similar data gaps as RHESSI and CGRO. HXRBS recorded a total of 12,772 flares above 
energies of > 20 keV during a life span of 9.8 years, so the average flare rate is < A >« 3.6 flares per day, 
which is about half that of BATSE and slightly less than RHESSI, which results as a combination of different 
instrument sensitivities, energy ranges, and flux thresholds. 

The waiting time distribution obtained from HXRBS/SMM is shown in Fig. 4 (top right), which has 
a similar range of At « 0.01 — 500 hrs as BATSE and RHESSI. We fit the same waiting time distribution 
(Eq. 18) with a powerlaw slope of p = 2 and find similar best-fit value for the average waiting times, i.e., 
At = 1.08 hrs, so the sensitivity and threshold is similar. Also the orbital effects modify the observed 
distributions exactly the same way. Thus we can interpret the distributions of waiting times for HXRBS in 
the same way as for BATSE and RHESSI, in terms of a nonstationary Poisson process with high intermittency. 

An earlier study on the waiting time distribution of HXRBS/SMM flares was published by Pearce et 
al. (1993), containing a subset of 8319 events during the first half of the mission, 1980-1985 (Fig. 4, top 
left panel). However, the published distribution contained only waiting times in a range of At = 1 — 60 
min (w 0.02 — 1.0 hrs), and was fitted with a powerlaw function with a slope of p = 0.75 ± 0.1 (Pearce et 
al. 1993; Fig. 5 therein). We fitted the same intermediate range of waiting times (At w 0.1 — 1.0 hrs) and 
found similar powerlaw slopes, i.e. p « 0.72 for HXRBS, p = 0.84 for BATSE, and p = 0.65 for RHESSI, 
so all distributions seem to have a consistent slope in this range. This partial powerlaw fits are indicated 
with a straight line in the range of At w 0.02 — 1.0 hrs in all six cases shown in Fig. 4. However, the waiting 
time distribution published by Pearce et al. (1993) extends only over 1.7 orders of magnitude, and thus does 
not reveal the entire distribution we obtained over about 5 orders of magnitude with HXRBS, BATSE, and 
RHESSI. If we try to fit the same waiting time distribution (Eq. 18) to the dataset of Pearce et al. (1993), 
we obtain a similar fit but cannot constrain the powerlaw slope for longer waiting times. 



- 9 - 



3.4. WATCH/GRAN AT Waiting Time Analysis 

A waiting time distribution of solar flares has earlier been published using WATCH/GRAN AT data 
(Crosby 1996), obtained from the Russian GRANAT satellite, launched on December 1, 1989. GRANAT 
has a highly eccentric orbit with a period of 96 hrs, a perigee of 2000 km, and an apogee of 200,000 km. Such 
a high orbit means that the spacecraft can observe the Sun uninterrupted without Earth occultation (i.e., 
no spacecraft night), which makes the waiting time distribution complete and free of data gaps. WATCH is 
the Wide Angle Telescope for Cosmic Hard X-Rays (Lund 1981) and has a sensitive detector area of 95 cm 2 
in the energy range of 10-180 keV. 

Waiting time statistics was gathered during four time epochs: Jan-Dec 1990, April-Dec 1991, Jan- April 
1992, and July 1992. Crosby (1996) obtained a waiting time distribution for n = 182 events during these 
epochs with waiting times in the range of At w 0.04 — 15 hrs. The distribution was found to be a powerlaw 
with a slope of p = 0.78 ± 0.13 in the range of At < 3 hrs, with an exponential fall-off in the range of 
At « 3 — 15 hrs. We reproduce the measured waiting time distribution of Crosby (1996, Fig. 5.16 therein) 
in Fig. 4 (middle left panel) and are able to reproduce the same powerlaw slope of p — 0.78 by fitting the 
same range of At « 0.1 — 2.0 hrs as fitted in Crosby (1996). We fitted also the model of the waiting time 
distribution for a nonstationary Poisson process (Eq. 18) and find a similar mean waiting time of Ato = 0.97 
hrs (Fig. 4), although there are no data for waiting times longer than > 15 hrs. Thus, the partial distribution 
measured by Crosby (1996) is fully consistent with the more complete datasets analyzed from HXRBS/SMM, 
BATSE/CGRO, and RHESSI. 

3.5. ISEE-3/ICE Waiting Time Analysis 

Another earlier study on waiting times of solar flare hard X-ray bursts was published by Wheatland 
ct al. (1998), which is of great interest here because it covers the largest range of waiting times analyzed 
previously and the data are not affected by periodic orbital datagaps. The International Sun-Earth/ Cometary 
Explorer (ISEE-3/ICE) spacecraft was inserted into a "halo" orbit about the libration point LI some 240 
Earth radii upstream between the Earth and Sun. This special orbit warrants uninterrupted observations of 
the Sun without orbital data gaps. ISEE-3 had a duty cycle of 70-90% during the first 5 yrs of observations, 
and was falling to 20-50% later on. ISEE-3 detected 6919 hard X-ray bursts during the first 8 years of its 
mission, starting in August 1978. 

Wheatland et al. (1998) used a flux selection of > 4 photons cm~ 2 s _1 to obtain a near-complete 
sampling of bursts, which reduced the dataset to n = 3574 events. The resulting waiting time distribution 
extends over At ~ 0.002 — 14 hrs and could not be fitted with a single nonstationary process, but rather 
exhibited an overabundance of short waiting times At < 0.2 hrs (Wheatland et al. 1998). We reproduce 
the measured waiting time distribution of Wheatland ct al. (1998; Fig. 2 therein) in Fig. 4 (bottom left 
panel) and fit the combined waiting time distribution of a nonstationary Poisson process (Eq. 18) and an 
exponential random distribution for the short waiting times (Eq. 1). We find a satisfactory fit with the two 
waiting time constants Ato = 0.03 hrs and At\ = 1.84 hrs. The primary component of Ato — 0.03 hrs (~ 2 
min) is much shorter than for any other dataset, which seems to correspond to a clustering of multiple hard 
X-ray bursts per flare, but is consistent with a stationary random process itself. The secondary component 
contains 45% of the hard X-ray bursts and has a waiting time scale of Ati = 1.84 hrs. This component 
seems to correspond to the generic waiting time distribution we found for the other data, within a factor of 
two. Since there are no waiting times longer than 20 hrs, the powerlaw section of the distribution is not well 
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constrained for longer waiting times, but seems to be consistent with the other datasets. 

4. Conclusions 

We revisited three previously published waiting time distributions of solar flares (Pearce et al. 1993; 
Crosby 1996; Wheatland et al. 1998) using data sets from HXRBS/SMM, WATCH/GRANAT, ISEE-3/ICE, 
and analyzed three additional datasets from HXRBS/SMM, BATSE/CGRO, and RHESSI. While the pre- 
viously published studies arrive at three different interpretations and conclusions, we are able to reconcile 
all datasets and the apparent discrepancies with a unified waiting time distribution that corresponds to a 
nonstationary Poisson process in the limit of high intermittency. Our conclusions are the following: 

1. Waiting time statistics gathered over a relatively small range of waiting times, e.g., < 2 decades as 
published by Pearce et al. (1993) or Crosby (1996), does not provide sufficient information to reveal 
the true functional form of the waiting time distribution. Fits to such partial distributions were found 
to be consistent with a powerlaw function with a relatively fiat slope of p « 0.75 in the range of 
At ~ 0.1 — 2.0 hrs, but this powerlaw slope does not extend to shorter or longer waiting times, and 
thus is not representative for the overall distribution of waiting times. 

2. Waiting times sampled over a large range of 5 — 6 decades all reveal a nonstationary waiting time 
distribution with a mean waiting time of At Q = 0.80 ± 0.14 hrs (averaged from WATCH/GRANAT, 
HXRBS/SMM, BATSE/CGRO, and RHESSI) and an approximate powerlaw slope of p « 2.0. This 
value of the powerlaw slope is consistent with a theoretical model of highly intermittently variations of 
the flare rate, in contrast to a more gradually changing flare rate that would produce a powerlaw slope 
of p w 3.0. Waiting time studies with solar flares observed in soft X-rays (GOES) reveal powerlaw slopes 
in the range from p w 1.4 during the solar minimum topw 3.2 during the solar maximum (Wheatland 
and Litvinenko 2002), which would according to our model correspond to rather gradual variation of 
the flare rate with few quiescent periods during the solar maximum, but a high intermittency with 
long quiescent intervals during the solar minimum. The flare rate essentially rarely drops to a low 
background level during the solar maximum, so long waiting times are avoided and the frequency 
distribution is steeper due to the lack of long waiting times. 

3. For the dataset of ISEE-3/ICE there is an overabundance of short waiting times with a mean of 
At = 0.03 hrs (2 min), which seems to be associated with the detection of clusters of multiple hard 
X-ray bursts per flare (Wheatland et al. 1998). 

4. The threshold used in the detection of flare events modifies the observed waiting time distribution 
in a predictable way. If a first sample of n\ flare events is detected with a low threshold and a 
second sample of events with a high threshold, the waiting time distributions of the form N(At) = 
(1/Atj)/(1 + At/Ati) 2 relate to each other with a reciprocal value for the mean waiting time, i.e., 
At2/Ati = (ni/ri2). This relationship could be confirmed for the RHESSI dataset for six different 
threshold levels. 

The main observational result of this study is that the waiting time distribution of solar flares is consis- 
tent with a nonstationary Poisson process, with a highly fluctuating variability of the flare rate during brief 
active periods, separated by quiescent periods with much lower flaring variability. This behavior is similar 
to the quiescent periods (soft state) and active periods (hard state) of pulses from accretion disk objects and 
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black- hole candidates (Mineshige et al. 1994). The fact that solar flares are consistent with a nonstationary 
Poissonian process does not contradict interpretations in terms of self-organized criticality (Bak et al. 1987, 
1988; Lu and Hamilton 1991; Charbonneau et al. 2000), except that the input rate or driver is intermittent. 
Alternative interpretations, such as intermittent turbulence, have also been proposed to drive solar flares 
(Boffetta et al. 1999; Lepreti et al. 2001), but they have not demonstrated a theoretical model that correctly 
predicts the waiting time distribution over a range of five orders of magnitude, as observed here with five 
different spacecraft. 
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Fig. 1. — One case of a stationary Poisson process (top) and four cases of nonstationary Poisson processes 
with two-step, linear- increasing, exponentially varying, and 5- function variations of the occurrence rate X(t). 
The time-dependent occurrence rates X(t) are shown on the left side, while the waiting time distributions 
are shown in the right-hand panels, in form of histograms sampled from Monte-Carlo simulations, as well as 
in form of the analytical solutions (given in Eqs. 10-18). Powerlaw fits N(At) oc At~ p are indicated with a 
dotted line and labeled with the slope p. 
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Fig. 2.— Top: Flare rate per day observed with RHESSI during 2002-2009, containing a total of 52,014 
events. Quiescent time intervals with At > 5 hrs are marked in the form of a "bar code" at the top of the 
panel. Bottom left: The frequency distribution of waiting times N(At) is shown for short time intervals 
At < 10 hrs, which shows peaks at subharmonics of the orbital period of « 1.6 hrs. Bottom right: The 
longer waiting time intervals At « 2 — 24 hrs can be fitted with a powerlaw function with a slope of p = 2.0. 
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Fig. 3. — Waiting time distributions for six different subsets of the data selected by thresholds, P > 
0,10,100,300,1000,3000 cts s" 1 . The full distribution (with no threshold) is fitted with a model of a 
nonstationary Poisson process (Eq. 18) with a powerlaw slope of = 2 (thick solid curve in top left panel). 
The same model functions arc predicted (using the scaling of Eq. 21) for the threshold-selected five subsets 
based on the number of detected events above the threshold. 
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Fig. 4. — Waiting time distributions of six different datasets: HXRBS/SMM (top left and right), 
WATCH/GRAN AT (middle left), ISEE-3/ICE (bottom left), BATSE/CGRO (middle right), and RHESSI 
(bottom right). The distribution of the observed waiting times are shown with histograms, fits of nonsta- 
tionary Poisson processes with dashed curves (with a powerlaw tail with slope of p = 2), and the best fit 
in the fitted range with thick solid curves. Powerlaw fits in the range of At w 0.1 — 2.0 hrs as fitted in the 
original publications (Pearce et al. 1993; Crosby 1996) are also shown (straight line and slope p). The excess 
of events with waiting times near the orbital period (t or ut ~ 1-6 hrs) is an artificial effect and is not included 



